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Abstract 

We show that the dynamics of an elastic solid embedded in a Minkowski space consist of a set of 
coupled equations describing a spin-1/2 field, ^, obeying Dirac's equation, a vector potential, A^, 
obeying Maxwell's equations and a metric, g^, which satisfies the Einstein field equations. The 
combined set of Dirac's, Maxwell's and the Einstein field equations all emerge from a simple elastic 
model in which the field variables ^f, A^ and are each identified as derived quantities from 
the field displacements of ordinary elasticity theory. By quantizing the elastic field displacements, 
a quantization of all of the derived fields are obtained even though they do not explicitly appear 
in the Lagrangian. We demonstrate the approach in a three dimensional setting where explicit 
solutions of the Dirac field in terms of fractional derivatives are obtained. A higher dimensional 
version of the theory would provide an alternate approach to theories of quantum gravity. 
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I. INTRODUCTION 



In constructing a quantum field theory, the usual prescription is to start with a known 
equation of motion, such as Dirac's equation, and "invent" a suitable Lagrangian that repro- 
duces the equation of motion when Lagrange's equations are applied. While this prescription 
has been successful it is not unique. In other words it is possible for two different Lagrangian's 
to lead to the same equation of motion. For example, section 12 — 2 of Reference[l] provides 
a good example of two different Lagrangian's that lead to the same equation of motion for 
the density variations in an acoustic field. 

In a theory of quantum gravity, this traditional approach would involve using a Lagrangian 
with an appropriate set of terms such that when Lagrange's equations are applied, the 
Einstein field equations are reproduced. The Lagrangian obtained in this manner explicitly 
contains the gravitational metric, as well as any other field variables that are coupled to it. 
For instance an attempt at merging gravity with QED would produce a Lagrangian that 
explicitly includes the Dirac Spinor field the electromagnetic vector potential and the 
gravitational metric g^ v . 

In this paper we demonstrate an alternate approach to quantum gravity based on a model 
of an elastic solid. In the this model, the only field variables that appear in the Lagrangian 
are the field displacements, Ui, that occur in elasticity theory. Using the methods of fractional 
calculus, we will show that the equations of motion of the system describe excitations that 
can be identified as massless, non-interacting, spin-1/2 particles obeying Dirac's equation. 
We then assume that one of our coordinates is periodic and use a dimensional reduction 
technique to reduce the dimensionality from three-dimensions to two-dimensions. 

When terms beyond the linear approximation are included, this dimensional reduction 
produces a new set of equations in which the spin field, is shown to interact via a 
vector potential and a metric with field variables g^ u . The compatibility equations of St. 
Venant are shown to reproduce Maxwell's equation for A^ and the Einstein field equations 
for g^ u . We quantize the field displacements using standard approaches and thereby produce 
a quantization of A^ and g^ u even though none of these quantities appears explicitly in 
the Lagrangian. We demonstrate the basic methods in a three-dimensional setting where 
exact expressions for the Dirac field can be obtained. 

When quantized, this theory provides a low dimensional version of a quantum description 
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electrodynamics coupled to gravity. If this procedure could be extended to higher dimensions 
it would provide an alternate approach to theories of quantum gravity. 

II. ELASTICITY THEORY 

The theory of elasticity is usually concerned with the infinitesimal deformations of an 
elastic body^^^. We assume that the material points of a body are continuous and can 
be assigned a unique label a. For a three-dimensional solid each point of the body may be 
labeled with three coordinate numbers a 1 with i — 1, 2, 3. 

If this three dimensional elastic body is placed in a large ambient three dimensional space 
then the material coordinates a 1 can be described by their positions in the 3-D fixed space 
coordinates x l with i = 1,2,3. We imagine that the solid is free to distort within the 
fixed ambient space described with coordinates x % . In this description the material points 
^(x 1 , x 2 , x 3 ) are functions of x. A deformation of the elastic body results in infinitesimal 
displacements of these material points. If before deformation, a material point a is located 
at fixed space coordinates x , x 02 , x 03 then after deformation it will be located at some other 
coordinate x 1 ,^ 2 ,^ 3 . The deformation of the medium is characterized at each point by the 
displacement vector 

u % = x % — x° l 

which measures the displacement of each point in the body after deformation. We will 
assume that our elastic solid is periodic in the coordinate a 3 and at various points in this 
paper we will Fourier transform the a 3 coordinate. 

It is one of the aims of this paper to take this model of an elastic medium and derive from 
it equations of motion that have the same form as Dirac's equation. In doing so we have 
to distinguish between the intrinsic coordinates of the medium which we will call "internal" 
coordinates and the fixed space coordinates which facilitates our derivation of the equations 
of motion. In the undeformed state we may take the external coordinates to coincide with 
the material coordinates a 1 = x 0t . The approach that we will use in this paper is to derive 
equations of motion using the fixed space coordinates and then translate this to the internal 
coordinates of our space. 
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A. Strain Tensor 



Let us assume that we have an elastic solid embedded in a three-dimensional Minkowski 
space with metric 

' -\ 
Vij = 10 
{ 1/ 

We first consider the effect of a deformation on the measurement of distance. After 
the elastic body is deformed, the distances between its points changes as measured with 
the fixed space coordinates. If two points which are very close together are separated by 
a radius vector dx° l before deformation, these same two points are separated by a vector 
dx % = dx° l + du l afterwards. The squared distance between the points before deformation 
is then ds 2 = —(dx 01 ) 2 + (dx 02 ) 2 + (dx 03 ) 2 . Since these coincide with the material points in 
the undeformed state, this can be written ds 2 = —(da 1 ) 2 + (da 2 ) 2 + (da 3 ) 2 = da l r]ijda J . 
The squared distance after deformation can be written 4 



ds' 2 = ^2 {dx 1 ) Vij (dx j ^ 

= ^2 (da 1 + du l ) rjij (dap + du j ) 



( du 1 \ ( du 3 \ 

= 5 r + ? a? ) mi r + ? s? J 

= ^2 riijda % da? + ^ r]ij——rda k da J + ^ r q i j—jda t da l + ^2 VijTrT~^rrda k da l 
■ ■ ■ ■ j (ja - j \ja --77 (ja (ja 

i J ijk i] l ijkl 

sr^( ( du k du k \ ^ du k du 1 \ 

= 5 r + ? v ik ^ + m ^ ) + ? Vkl d^^ ) da da 

= y2 (.Vij + 2%) da 1 da 3 

ik 

where is 

1^/ <9w fe <9u fc ^ <9-u fc du 1 \ 

e - = 2 ? (** + ^ + E J • ( 2 ) 

and the presence of the matrix ^ simply reflects the fact that we are assuming our solid is 
embedded in a Minkowski space with a pseudo-Euclidean metric. 
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The quantity is known as the strain tensor. It is fundamental in the theory of elasticity 
In the above derivation, the material or internal coordinates were treated as functions of the 
fixed space coordinates. As is well known in elasticity theory we could just as well treat 
the fixed space coordinates as functions of the material coordinates. In this case, the strain 
tensor has the form 3 



These two different approaches to the strain tensor are known in elasticity theory as the 
Lagrangian and Eulerian perspectives. In this work we will derive the equations of motion 
using the fixed space coordinates which simplifies the derivation and we will translate the 
result, when necessary, to the internal coordinates. 

In most treatments of elasticity it is assumed that the displacements u l as well as their 
derivatives are infinitesimal so the last term in Equation (j2j) is dropped. In this work, we 
will treat the strain components as small but finite. We will then examine the structure of 
the equations of motion when the higher order terms are treated as small perturbations on 
the infinitesimal strain results. 

B. Metric Tensor 

The quantity 



is the metric for our system and determines the distance between any two points. One in- 
teresting aspect of the elasticity theory approach is that it provides a natural metric on the 
system in terms of the strain components expressed entirely in terms of the internal coordi- 
nates of the elastic body. This means that at any point in space the distance measurement 
can be made without reference to the fixed space coordinates. In other words if you were an 
ant living in this elastic medium, Equation (J3J) would be the metric that you would use. 

Even though the metric in Equation (j3J) does not have the Euclidean form, the space in 
which we are working is still intrinsically flat. The metric that we derived is due simply to 
a coordinate transformation and so cannot describe the curved space of general relativity. 




(3) 





Vij 
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That this metric is simply the result of a coordinate transformation from the Minkowski 
metric can be seen by writing the metric in the forrn^ 

( -1 q\ 



( dx 1 


dx 2 


dx 3 \ 
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da 1 


dx 1 


dx 2 


dx 3 


da 2 


da? 
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dx 1 
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dx 3 


V da 3 
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da 3 1 
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/ dx 1 


dx 1 


dx 1 \ 
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da 3 


dx 2 


dx 2 


dx 2 


da 1 


da 2 


da 3 
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da 3 I 



where 



dx^ 
da v 



J T T}J 



du^ 
+ da^- 



and J is the Jacobian of the transformation. Later in section IVIl however we will use a 
dimensional reduction technique borrowed from Kaluza-Klein theories to reduce the three- 
dimensional flat space to a two-dimensional curved space. We will show that the metric for 
the Fourier modes of this two dimensional system is not a simple coordinate transformation. 

The inverse metric, which is written with upper indices as g lk , can be obtained by explicitly 
inverting Equation (j3J) or we can write {g lk ) = (J^ 1 )i]^ 1 (J^ 1 ) T where 

/ da 1 da 1 da 1 \ 



J- 
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dx 3 
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dx 3 
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da 3 



(5) 



V dx 1 dx 2 dx 3 / 



This yields for the inverse metric 
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fJtU 



du 11 
ua dx^ 



du v ^ du^ du v 

a/3 



(6) 



Equation ([HD shows that we can write the inverse matrix directly in terms of derivatives of 
Ui with respect to the fixed space coordinates. This form of the inverse metric will be useful 
in later sections. 



1. Internal vs. External Coordinates and Summation Convention 



The change in form of the metric between that given in Equation ([T]) and Equation (j4j) is 
due simply to a change in coordinates between the fixed space coordinates and the material 
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coordinates. In this regard the transformation is similar to changing from Cartesian to 
spherical coordinates. This change is useful because it allows us to derive equations in the 
fixed space coordinates where the calculations are simplified, and then when necessary we 
can switch to the internal coordinates using u l = x l — a 1 . 

We would like to be able to use the notation that a raised index on a variable indicates a 
contraction with the metric tensor and that a raised index and a lower index with the same 
label implies a summation (ie the Einstein summation convention). We have to be careful, 
however, to point out which set of coordinates, and hence which metric we are using, so we 
will be explicit in each section as to which coordinate system the raised indices refer to. For 
instance we can write Equation ([2]) more compactly as 

1 / dui duj dui du l \ , 

where uiu 1 = J2k VikUiu k so that upper/lower indices indicate contraction with the fixed space 
metric, Equation (JTJ. 



C. St Venant's Equations of Compatibility 

In Section flHDI) we will derive the equations of motion of the elastic solid using the 
Lagrangian formalism. There are, however, additional constraints that an elastic solid must 
also satisfy. These constraints are called the St Venant equations of compatibility in classical 
elasticity theory™. The usual description of these compatibility equations is that they are 
integrability conditions or, a restriction on the strain components such that they can be 
considered partial derivatives of a function u as displayed in Equation (j2J). In other words if 
6ij is a function that is composed of the partial derivatives of u then it has to satisfy certain 
conditions and these are the compatibility equations. 

However, from a geometric standpoint these equations are simply a re-statement of the 
fact that space is flat£. In other words the compatibility equations are equivalent tc 6 

r - 1 I d2gav d2g ^ d 2 gp u 

aPflu - 2 y da p^ da(iQ y da a aU daa ^ 

where R a pn V is the Riemann Curvature tensor and T p a/3 are the Christoffel symbols given 

by£ 

f p x P ( ^9\a dgxp _ dg aP \ 
af3 9 { dx? dx a dx x ) ' 
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In the above equations an upper/lower index implies a contraction with the metric in the 
internal coordinates. 

One of the interesting aspects of an elastic solid is that this setting gives you for "free" 
an explicit expression for the metric, Equation and a statement about the curvature of 
space, Equation (jSJ). We get these equations even though, as we will see shortly, the metric 
is not a dynamical variable appearing in the Lagrangian. 



D. Equation of Motion 

In the following we will use the notation 



du L 



and therefore the strain tensor is 



]- (u^ + u UfM + u ¥ u l v ) . (9) 



e M „ 2 

and all contractions are with the fixed space metric, Equation ([!]). 

We work in the fixed space coordinates and take the negative of the strain energy as the 
lagrangian density of our system. This approach leads to the usual equations of equilibrium 
in elasticity theory^. The strain energy is quadratic in the strain tensor e^ u and therefore 
the Lagrangian can be written 



L = ~ E C^e^e 



ap 



puap 



The quantities C ,ll ' ap are known as the elastic stiffness constants of the material. For an 
isotropic space most of the coefficients are zero and in fact there are only two independent 
elastic constants in a three-dimensional isotropic space. The lagrangian density then reduces 
to 

L = (A + 2/i) [e^ + e\ 2 + e^] + 2A [e 22 e 3 3 - ene 2 2 - £11633] + 4/i [e^ - e\ 2 - e? 3 ] (10) 

where A and /x are known as Lame constants^. 

We first derive the equations of motion of the system in the approximation where the strain 
components Uij are infinitesimal. In the infinitesimal strain approximation, the quadratic 
terms in Equation ([9]) are dropped giving 
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The usual Lagrange equations, 

V — (-^-\ — 

V dx u \du pv J du p 

apply with each component of the displacement vector, u p , treated as an independent field 
variable. 

Using the above form of the Lagrangian one can write 

dL 



du pv 



2X(arj pv ) + 4/xe 



pv 



where the divergence of the displacement field is o = (— u n + ^22 + ^33)- In classical elasticity 
theory a is known as the dilatation and physically represents the fractional change in density 
of a medium due to a deformation. 

We now have three field equations (one for each value of p), 

(2A + 2 y u) T f-a + 2pV 2 K) = (11) 



dx u du pv dx p 

where V 2 = — ^ + ^ + J^-. Applying the operator d/dx p to Equation (TlTi) yields the 
wave equation 

/ d 2 d 2 d 2 \ 

(-^ + S5 + &l)^° (12) 

Equation ( fl2l) shows that the classical dilatation in the medium obeys the wave equation. 
In Section ( 1IIII) we will demonstrate a new method for reducing the wave equation (fl2l) to 
Dirac's equation and compare this method to the traditional Dirac reduction. But first we 
turn our attention to quantizing the field displacements Ui in this elastic model. 



E. Quantization 

The Lagrangian density of the system is given by Equation ( flOl) . The coordinate X\ plays 
the role of time in this three-dimensional space so the canonical momenta associated with 
the field variable u p are given by 

p du 



pi 



which gives 



P 1 = -2A(e n - e 22 - e 33 ) - 4/xen (13) 
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P 2 = 4//612 

Pa = 4/1613- 

The Hamiltonian density is denned as 

n = £i>?-L 

p 

= -P\Uu + P 2 «2i + -P3M31 - £ 
and the total Hamiltonian is the integral over all space of the Hamiltonian density 

H = I d 3 xH 



Inverting Equation (fl3l) allows us to replace the variables u p % with P p in the result. This 
gives 



n _ 3 A 2 Pj + Pj A (u 22 + ^33) P 1 

U - 4(A + 2^ + - ^ AT2^ P2Ul2 " P3Ul3 

+/•* ( M 23 + M32) 7- — h 2 \U 2 2U 33 

+(A + 2/i)(^ 2 + M y (14) 

We now Fourier transform the field variables making the assumption that one of our coor- 
dinates, a 3 is compact with the topology of a circle. Therefore, when we Fourier transform 
the field variables, the #3 component is associated with a discrete spectrum while the other 
two coordinates are continuous. Writing out the coordinate dependencies explicitly we have, 



P p (x 1 ,x 2 ,x 3 ) = J2 dqx dq 2 P p ^ iq ' x 

93 J J 

u p (x 1 ,x 2 ,x 3 ) = 53 J dqi J dq 2 u p ^ tq ' x . 



93 

The Fourier transform results in terms in the Hamiltonian that mix field variables associated 
with q and —q. For instance, the contribution to the total Hamiltonian from the Pj 2 term in 
Equation (TT41 becomes 



j d 3 xP? = J d 3 x (y, J dqi J 


fdq 2 Y.J 


dq[ j 




<?3 





E / dqi [ dq 2 P 1;? P lr 
93 J J 
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The total Hamiltonian then becomes 



i 

where H q is written symmetrically in q and — q as 



q " A + 2/i A + 2/x A + 2/i A + 2/i 1 J 

+ 2(A + 2/i) 2(A + 2/i) + fe^^ + 2(A + 2/i) 
ig3-Pi _„M3 „A 2 

+^293^2,3^3 ,- g A 2(A + 2/i') + q2q ^ U2 -1 U ^ X + % U 3~q U 3,q X 

. 3Pl ,-qP\,q P2,-qP2,q . Pz,-qP%,q . 1. D 

+ 4(AT^) + "^T" + ~V~ + 2^^- (16) 

+fiqlu 2 ,-qU 2) g + //?2?3«2,g"U3,-g + jU?2?3«2,-9«3,g + ^2 U 3,- g «3,g + 2/ig|u 3 -9«3,g (17) 

Since terms in the Hamiltonian with different values of g are not mixed, the function # 5 in 
Equation ( TT5l) can be solved independently for each g. if g is a bilinear function in the vari- 
ables w^ig and Pi,± g and can be diagonalized exactly, using the methods of Biougliobov^^. 



1. Exact Diagonalization 

The methods used in diagonalizing the Hamiltonian in Equation (|15[) are summarized in 
the references' 1 ^. The idea is to rewrite the Hamiltonian in terms of a set of creation and 
annihilation operators, bi >q and b iq such that the Hamiltonian has the form 



l,q 



H q = J2 ^i,q h lq h 
i 

and the operators satisfy the commutation relations 

lKqi b ],qi] = il5 i,jSq,q'- 

The details of this procedure are included in Appendix (jAl) . Of particular interest are the 
energy eigenvalues of the modes. There are three distinct positive energies for the states b{. 
They are 

1 
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1 



4/x(A + fj) 



(19) 



4 



^ y / /i(A + /i)(A + 2 yU ) 



E; 



1 



4// (A + //) 



(20) 



4 



With the operators bi calculated, the field variables Ui can be written as a linear combi- 
nation of creation and annihilation operators as 



where the c^- are coefficients given in the appendix and q\ is the energy of a given mode. 
These eigenstates are the linear approximation obtained from keeping the lowest terms in 
the strain components Uij in Equation (TTU]) . The higher order terms that were left out of the 
Lagrangian can be incorporated by treating them as perturbations. In other words we can 
use standard perturbation theory to find new strain components that are nonlinear in the 
creation and annihilation operators. These field components can then serve as the basis for 
a theory of finite strain as we do in the next section. 

Strictly speaking the field displacements as expressed in Equation f[2"Tj) are not energy 
eigenstates since the fej have different energies. We will mainly be concerned however, with 
a low energy approximation of the spectrum of this elastic solid. For positive values of A 
and /i and A > (6 + 4a/3^ /z, the energies E 3 q can be arbitrarily small compared to E X q 
and E 2t q. For instance with A = 13/i and /i = .1, the energies Ei q and E 2tq are more than 
20 times greater than E% yq . This suggests that in a low energy theory only the excitation 
corresponding to energies E 3 ^ q will be present for suitably defined Lame constants. We will 
not investigate the mechanical properties of such a solid but merely point out that in such a 
theory at low energies, both the field displacements Ui and the dilatation, V • u, are energy 



eigenstates. We also note that each of the energies is proportional to J q 2 + q$ with q 3 taking 
on discrete values. So one would expect in the lowest energy approximation that only the 
modes with q 3 = will be present and at slightly higher energies the mode with q 3 = 1 will 
be present. This low energy approximation will be exploited in later sections. 

One of the things that we have gained from this formalism is the ability to calculate any 
quantity, that depends on the field displacements, quantum mechanically. For instance we 
can now calculate the metric given in Equation ([6]) using the form of the field decomposition 
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(21) 



q=l i=l 
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given in Equation (}2~Tj) even though the metric itself is not a dynamical variable appearing 
in the Lagrangian. 

In the finite strain theory treated in Section (JVj we will need to take Fourier transforms of 
field variables in the internal coordinates rather than the fixed space coordinates. Since we 
will be keeping the nonlinear terms in all of our equations, then for consistency we assume 
that the field variables in Equation (I2T!) have been properly treated to the same order in 
perturbation theory. We will not explicitly calculate the field variables in perturbation 
theory rather we will focus on the form of the field equations when terms beyond the linear 
approximation are kept. 

We will now give a new derivation of Dirac's equation as the equation of motion of the 
elastic solid. 

III. DERIVATION OF DIRAC'S EQUATION OF MOTION 
A. Cartan's Spinors 

The concept of Spinors was introduced by Eli Cartan in 1913^. In Cartan's original 
formulation spinors were motivated by studying isotropic vectors which are vectors of zero 
length. In three dimensional Minkowski space the equation of an isotropic vector is 



for generally complex quantities X{. A closed form solution to this equation is realized as 



The two component object £ = (£o,£i) has the rotational properties of a spinor 12 and any 
equation of the form (1221) has a spinor solution. 

In the following we use the notation = d/dx^ and the wave equation is written 



(22) 



xi = Co + f ?, x 2 = £o - £i> and ^3 = -2f f i 



(23) 



where the two quantities & are then 




{-dl + d\ + dj) = 0. 
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This equation can be viewed as an isotropic vector in the following way. The components of 
the vector are the partial derivative operators d/dx^ acting on the quantity <p. As long as 
the partial derivatives are restricted to acting on the scalar field it has a spinor solution 
given by 

^ = 2 (afci + dx*) (24) 

and 

^ = 2 [dx* ~ dx^j (25) 
where the "hat" notation indicates that the quantities £ are operators. Let us now introduce 
the variables 

z = x\ + x 2 and Z\ — X\ — x 2 
Equations (1241) and (|25p are now 

so 



and 



dz° 



dz 1 

These are equations of fractional derivatives of order 1/2 (also called semiderivatives) denoted 
<fo = D l J? and fi = D 1 ^ 2 . Fractional derivatives have the property that 13 

D l ' 2 D 1 ' 2 = — 

z z dz 

and various methods exist for writing closed form solutions for these operators^ 1 ^. The 
exact form for these fractional derivatives however, is not important here. The important 
thing to note is that a solution to the wave equation can be written in terms of spinors which 
are fractional derivatives. 

One of the interesting properties of suitably defined fractional derivatives (for instance 
the Weyl fractional derivative) that will be exploited in later sections is their action on the 
exponential function. While the derivative of an exponential is given by 

d 

e ax = ae ax 

ox 

the semiderivative of the exponential function is given by 

D l J 2 e ax = y^e ax (26) 
This will prove useful later when we Fourier transform the equations of motion. 

14 



B. Matrix Form 



It can be readily verified that our spinors satisfy the following equations 





dx 3 



+ 6 



d d 
+ 



dx 1 dx 2 



a 



d 



and in matrix form 



dx 2 dx 1 

( d d 

dx 3 

d d 



+ 



d 

dx 3 



d \ 



The matrix 



V dx 2 dx 1 
( d 

X 



dx 1 dx 2 

d_ 

dx 3 > 





dx 3 

d 



+ 



= 
= 

(t \ 

d \ 



(27) 



dx 1 dx 2 
d d 



V dx 2 dx 1 dx 3 ) 

is equal to the dot product of the vector <9 M = d/dx^ 1 with the pauli spin matrices 



X 



d 



■7 2 + 



d 



dx 1 dx 2 dx 3 

where 

-l o i V 10 / y o — l 

are proportional to the Pauli matrices and satisfy the anticommutation relations 



I 



\ 



{ 7 ^} = 2/7/ 



(28) 



where I is the identity matrix. 
Equation (|27|) can be written 

]T^ 7 ^> = 0. (29) 

This equation has the form of Dirac's equation in three-dimensions for a noninteracting, 
massless, spin- 1/2 field, 



C. Relation to the Dirac Decomposition 

The fact that the wave equation and Dirac's equation are related is not new. However 
the decomposition used here is not the same as that used by Dirac. The usual method of 
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connecting the second order wave equation to the first order Dirac equation is to operate on 
Equation fl29|) from the left with J2l=i l v d v giving 



3 



= E l u ^d v d^{x) 



//,f=l 



= (-dl + dl + dl) *(x) 



(30) 



where ^ = (a%, a^) is a two component spinor and Equation (128]) has been used in the last 
step. 

This shows that Dirac's equation does in fact imply the wave equation. The important 
thing to note about Equation (1301) however, is that the three-dimensional Dirac's equation 
implies not one wave equation but two in the sense that each component of the spinor \I/ 
satisfies this equation. Explicitly stated, Equation fl30|) reads 



for the independent scalars ai,a 2 . 

Conversely, if one starts with the wave equation and tries to recover Dirac's equation, it 
is necessary to start with two independent scalars each independently satisfying the wave 
equation. In other words, using the usual methods, it is not possible to take a single scalar 
field that satisfies the wave equation and recover Dirac's equation for a two component 
spinor. 

What has been demonstrated in the preceding sections is that starting with only one 
scalar quantity satisfying the wave equation, Dirac's equation for a two component spinor 
may be derived. Furthermore any medium (such as an elastic solid) that has a single scalar 
that satisfies the wave equation must have a spinor that satisfies Dirac's equation and such 
a derivation necessitates the use of fractional derivatives. 

IV. DIMENSIONAL REDUCTION IN INFINITESIMAL STRAIN 

In this section we take a closer look at the equation of motion, Equation (1291) . and the 
spinor (represented as a fractional derivative) when the field displacements are Fourier trans- 
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formed. In the infinitesimal theory of elasticity, all terms in mj beyond the linear term are 
dropped. In the infinitesimal theory therefore, no distinction is made between transforming 
the coordinates a; and Xj. Later when we assume small but finite strain components, we will 
need to distinguish these coordinates. 

The Dirac field in Equation (|2T?|) is given explicitly by 

/ 



i(j)(xx, x 2 ,x 3 ) 



Dll 2 
Dl( 2 



0(xi,x 2 ,x 3 ) (31) 



with z = X\ + x 2 and z\ = X\ — x 2 . Because the semiderivatives -D^ 2 Zl are independent of 
X3 we can bring e lXs through the operator £ when we Fourier transform <p. 

Transforming the dilatation first in the periodic coordinate (73, Equation fl29|) becomes 



= J2Yd,i<p 
= E e * 93X3 ( E 7^ + 7^3) ^ 

= E ^ (e TVfy + iq, j |0 93 (32) 

where we used = J2q 3 (pq3 emX3 , an d 7 3 7 3 = 1- 

Equation ( 132]) is equal to zero only if the coefficients of e mx:i are zero for each value of q 3 
giving 

EH'^-ga W = (33) 

where 7 //x = 7 3 7^ and satisfies the conditions for a two dimensional metric {7'^, r y' u } = llvj^ 
with 

mi/ 



/1 



(34) 



-1 

Equation (1331) shows that the fourier modes of the elastic solid obey a two dimensional version 
of Dirac's equation for spin-1/2 particles with a mass q 3 . The two continuous variables left 
in the problem are X\ and x 2 with X\ playing the role of time. 
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Let us examine the form of the spinor further by Fourier transforming the two continuous 
coordinates x\ and x 2 . Using Equation fl26l) we further transform the spinor as 



£(j} q3 (x 1 ,x 2 ) 



D l J 2 

° | ^3(^1^2) 
D\' 2 




rtl/2 i(q 1 x 1 +q 2 x 2 ) 
Jjl/2 e 'i(qixi+q2X2) 



J q\ mm 



£)l/2g|izo(i?l+<?2)+§ iz\{q\-qi) 

Jjl/2 e \iz Q (qi+q2) + \tz 1 ( y qi-q2) 
2*1 



J Qimmi 



\fi<p. 



91.92,93 



(35) 



d qi dq 2 e^i^+^) 

When the fields fej are viewed as a time dependent quantity (ie the Heisenberg picture 
in quantum mechanics), the wavevector q\ is equal to the energy of the q th mode^£. This 
allows us to write the column vector in Equation (1351) as 

/ 

u{q) = 



VE + q 2 



y y/E - q 2 

Compare this formula to the expression for the four dimensional Dirac spinor in a definite 
state of helictyi£ 



u(p) 



( 


^0^ 


\ 


VE + P 






I 1 J 




VE-P 


f 0^ 




V 


11 / 


/ 



The dimensional reduction in the infinitesimal theory of elasticity has produced a two 
dimensional version of Dirac's equation for a particle with a bare mass q% and a spinor that 
has a consistent form to the known 4 dimensional version. 

Incidentally, with the identification q 3 as a mass term, the energy eigenstates given in 
Equation f|T8|) are seen to have the relativistic form E q ~ J ' q\ + m 2 . 

The only problem with the interpretation as relativistic particles is the question of quan- 
tum statistics. Relativistic spin- 1/2 particles obey anticommutation relations while the bo- 
son operators hi that we have defined obey commutation relations. This situation could be 
"rescued" by defining new operators 

Cfc = k b k 



where the 6 k are complex Grassman numbers. Grassman numbers satisfy 

Ok@k' = —Ok'Ok, (9k) 2 = 

this would imply 

{cfc,c fc /} = {0 k b k ,8 k >b k >} 

= 9 k 9 k 'b k bk> + k >9 k b k ibk 
= OkQk'bkbk' — d k O k 'b k 'b k 
= Q k Q k i [b k , by] 
= 

Similarly we have, 

{ c L c l} = °> and { c > c l} = i>QkQi>h,k' 

With this definition, the fields c k satisfy appropriate quantum statistics and still obeys 
Dirac's equation. We do not wish to dwell on this admittedly ad-hoc procedure for getting 
fermion statistics into this theory. Rather we wish to focus on the form of the equations of 
motion that are derived and demonstrate that the all ingredients for a theory of quantum 
gravity are present in this model. 

In this section, we have demonstrated that the equation of motion of this model of an 
elastic solid, in the infinitesimal strain approximation has as its equation of motion a two 
dimensional version of Dirac's equation for spin-1/2 particles. However, these spin-1/2 par- 
ticles described in Equation (1331) do not interact with each other or the outside world. In the 
next section we will treat the dimensional reduction problem again in the context of the finite 
theory of strain. When terms beyond the linear approximation are kept, we show that the 
equations obtained describe particles that interact gravitationally and electromagnetically. 

V. FINITE STRAIN 

A. Internal Coordinates 

In this section we will need to Fourier transform our field variables in 03 and therefore 
need to translate the equations of motion from the fixed space coordinates to the internal 
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coordinates. For clarity and to adopt a more consistent convention, in the remainder of this 
text we change notation and write the internal coordinates not as a 1 but as x n and the fixed 
space coordinates will continue to be unprimed and denoted x l . Now using u % = x l — x' % we 
can write 



d ^ dx' j d 



dx % j dx l dx'j 



E 



' dxi du? \ d 



j 



E * 



dx % dx % J dx'j 
dv? \ d 



j 



lJ dx i / dx'j 



3 dv? 3 

E aZi oZm ft®} 



dx H ^ dx i dx'j 

3 



Equation fl36l) relates derivatives in the fixed space coordinates x l to derivatives in the 
material coordinates x H . We can now re-write the three-dimensional Dirac's equation as 



3 3 / du v \ 

E-W = Erk-E^ * 



= E = (37) 

11=1 

where ^ = £0, 9' = d/dx'^ and 7'^ is given by 

y/. = 7 m _ £ <7 <*. (38) 

The 7'^ are simply the gamma matrices expressed in the primed coordinate system. The 
anticommutator of these matrices is 

{7' M ,f} = {7 M -E<7 a ,7 , '-E«V} 

a /3 

= {7 M , 7"} - E ^{7 M , 7^} " E <{7 a , 7*} + E <^{7°, 1 P } 

[3 a a/3 



21 ir - E - E « V" + E E «» c 



a a fj 



Comparison with Equation ([6]) shows that 

{ 7 "\7"} = 2Ig^ 
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(39) 



These gamma matrices have the form of the usual dirac's matrices in a curved spaced. 
To further develop the form of Equation fl37|) we have to transform the spinor properties 
of £. Much like a normal vector, the components of a spinor are altered under a change of 
coordinates and as currently written £ is a spinor with respect to the Xj coordinates not the 
x\ coordinates. To transform its spinor properties we assume (similar to Brill and Wheler— ) 
a real similarity transformation and write \I/ = S^' where S is a transformation that takes 
the spinor in x^ to a spinor in x 1 . 

We then have 

Equation fl37|) then becomes 

= i^Sd'^' + (d'^S)*'} 
= ^S[d'^' + S~\d^S)^'] 
= S-WSldft' + S-^d'^S)^'} 
Using (d'^S^^S = — S 1-1 ^ S 1 ). This can finally be written 

7T;-rg*' = o (40) 

where T M = (d'^S^^S and 7 M = S~ lr -f ,fl S. The new gamma matrices 7^ still satisfy the 
appropriate anticommutation condition 

{r,l u } = 2Ig^ (41) 

What we have done in the above manipulations is re-write Equation (129]) completely in 
terms of the internal coordinates of the elastic solid. Equation (l4"Uj) has the same physical 
content as Equation ( |29|) . it is just expressed in different coordinates. 

Notice however that Equation fj4"0l) now superficially has the form of the Einstein-Dirac 
equation in three-dimensions for a massless, noninteracting, spin-1/2 particle^. The quantity 
9^— T M is the covariant derivative for an object with spin. In order to make this identification, 
the field must satisfy the additional equation^^ 

^7 + 7^-^ + 7^ = 

where Ya v is the usual Christoffel symbol. That this equation holds is shown in Appendix flBj) . 

Just as was done in the infinitesimal strain approach we will use a dimensional reduction 
to transform Equation (|40j) . In this case however, the gamma matrices couple the metric 
into the problem and this metric will also need to be reduced from three to two dimensions. 
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VI. DIMENSIONAL REDUCTION IN FINITE STRAIN 



The dimensional reduction method that we use is borrowed from Kaluza Klein theory. It 
consists of two parts. First we need to remove the dependence of the field variables on the 
coordinate x' 3 . This is accomplished by Fourier transforming the field variables as before. 
The second part consists of reducing the 3x3 matrix, g^ V) in three dimensions to a 2 x 2 
matrix suitable for a two dimensional space. It is this reduction of the metric that will 
introduce the electromagnetic vector potential into our equations. 

Let us denote the 3x3, three dimensional metric as g^ u and the 2x2, two dimensional 
metric as g^ v . We now write the Kaluza Klein ansatz^22 



I g a/3 - $ 2 A n -$ 2 A„ \ 



9 a/3 



I g a ? -A a \ 



(42) 



where the vector A^ is the electromagnetic vector potential, $ is a scalar and the upper/lower 
indices indicate a contraction with g^ v . 

With this ansatz for the metric, it can be shown that the flat space condition R a p = in 
three dimensions, reduces to the Einstein Field equations in curved two-dimensional space 
and Maxwell's equations for the field A^. The details of this derivation are given by Liu and 
Wesson 2 ^ where the reduction is performed in a five dimensional setting. This is a standard 
dimensional reduction used in Kaluza Klein theories and has been shown to produce a 
consistent treatment of four- dimensional gravity and electromagnetism from five-dimensional 
flat space. Therefore, we will not repeat this derivation here only noting that there is nothing 
in the derivation that is particular to five dimensions. The exact same treatment works in a 
dimensional reduction from three to two dimensions so that we will freely quote the results of 
that work which gives the following equations for the metric and the electromagnetic vector 
potential 



R a p = --$ 2 F A ^F /3 A + $- 1 $ Q;/3 



where F a p = A^ a — A a -p and we have used the notation where a comma indicates an ordinary 
derivative and a semicolon indicates covariant differentiation. 
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The second of these equations, may be written in the more traditional form 20 
where the quantities and T S Q/3 are effective energy momentum tensors given by 

rpa/3 _ ^-1 ^a;P _ g a P<§>^ 

So what we see is that by using the dimensional reduction technique of Kaluza Klein theory 
we automatically obtained the Einstein field equations and Maxwell's equations. What we 
now show is that the same dimensional reduction technique not only produces these equations 
but also changes the non-interacting Einstein-Dirac equation in three dimensions, into an 
interacting theory for a massive spin-1/2 particle in two dimensions. 



A. Dimensional Reduction of the Dirac Equation 

We begin by rewriting Equation (|40p as 

E r K - rvi + f & - r s ] ^ = o (44) 

We now rewrite the matrix 7 3 . From Equations ()42p and (j4ip we can write 

7 3 ^ 3 = g 33 = A^A^ - $~ 2 

This allows us to write 

T^Ef^ + Ti*" 1 (45) 

o=l 

where the matrix 7j_ satisfies 

{7 1 ,7±} = 0, {7 2 , 7± } = 0, 7 i = l 

An explicit expression for 7j_ is given in Appendix (jC|) . It can be seen that this form for 7 3 

correctly gives 7 3 7 3 = g 33 . We can now write Equation (1441) as 

2 2 

= E7 M [^-r,]vi/' + (^7^ a + 7^- 1 )K-r3]^ 
= E + A M - r 3 ) - r M ]<f + ^[d', - r 3 ]*' 
= E 7'"K + ^ - r 3 ) - r M ]V + s^ft - r 3 ]*' (46) 
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where ^"^ = 7_i_7 At . This equation has the same form as the classical Einstein-Dirac-Maxwell 
equation^^i 1 ^ 1 ^ 3 -. The metric enters in the equation through the j" 11 matrices which satisfy 
{y/^y^j = 2Ig^ v and the electromagnetic vector potential enters in the same way as the 
minimal coupling prescription <9 M — > d^ + ie/cA^. Additionally, a mass term has been created 
in Equation (HD with m = <I> _1 [c% - T 3 \. 

The field variables in Equation (l4"6{) still are dependent on x 3 . To complete the derivation 
we need to remove this dependence by Fourier Transforming our field variables. We first 
write the spinor field 

i 

this gives 

E e%qx ' 3 1 E ^""K + - r 3 ) - rg*' + ®- l [iq -r 3 ]A = o 

<? U =1 J 

Now we write 
Which gives 

q,k,k',k" [fj,=l 

This equation is true only if the coefficients of the exponential are independently true. Let 
m = q + k + k' + k" then 

E { E 7H<^A,cA",o + A^ k (i(m -k-k' - k")5 k » fi - r 3 , fc ») - ^ ^,kh" flWm-k-k' -k" 

k,k',k" [fJ.=l 

+5 kll ^- k }[i{m -k-k 1 - k")5 kfi - r 3>k }V m _ k _ k ,_ k „} = 0. (47) 

This is a series of equations, one for each distinct value of m. The dynamics contained 
in Equation (14?]) describes an infinite series of spin particles ^ m ^ k ^ k i^ k n interacting via a 
infinite series of vector potentials A^ k . Up until now we have kept all terms in this series of 
equations. We will now examine Equation (T4T|) in the low energy approximation. 
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B. Spectrum of Lowest modes 



Each of the quantities in Equation (j47j) is a function of the field displacements itj. We 
showed in previous sections that the energy of the Fourier modes of field displacements 
increases with q 3 . We therefore expect that in a system where the energy is arbitrarily low, 
not all of the modes in Equation ([17)1 will be excited. At the lowest energies only the mode 
u q=Q will be excited, as the energy of the system increases modes u q= ±\ becomes excited 
and so on. As this energy is increased Equation (T4"7|) rapidly becomes more complex but at 
sufficiently low energies its form is quite simple. 

We illustrate this by considering a theory in which only the lowest modes are present. 
Strictly speaking, our Hamiltonian formalism in Section (III E[) produced energy eigenstates 
for modes u q in a basis e iqx while we want to make a statement about the energy of the 
modes,Mq/ in the basis e tq ' x \ with x H = x l — u 1 . In Appendix ((D|) we show that the two are 
related and that if the only mode present in the system is q = in the first basis then, the 
only mode present is q' = in the second. 

Let us now consider the case where there is insufficient energy to excite the mode u q =±i 
and only u q=0 is excited. This would represent the lowest energy possible in our system. 
In this case the wavevectors in Equation (1471) must equal to m = k = k' = k" = and 
Equation (T4"7j) reduces to 

£ 7oT; - A,, r 3 , - r^o]^ - $0^3,0^0 = 0. (48) 
From Equation (jH]) we also have 

g*» = rT + il \-r] ua ^- -V" a ^ + y V a ^^4) , (H,v= 1, 2, 3) 
y I dx a 1 dx a £i ' dx a dx? J v ' ' ; 



which via Equation ( 1421) . provides an explicit expression for the two dimensional metric and 
the electromagnetic vector potential 

r = r + g ^ ^ _ ^ + g I , 2) (49) 

and 
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and where 

( -I q\ 



I • (51) 

V o o i; 

Note that the 2D metric in Equation (H9|) is no longer due to a simple coordinate change. 
There is no coordinate transformation (or any other transformation) that involves only the 
coordinates x\ and x' 2 that will remove the Fourier transforms in this equation and globally 
create the form given in Equation ([T]). 

Let us summarize what we have done in this section. We began with an equation for 
Dirac's equation in three-dimensions, Equation (|40j) . which described a spin-1/2 particle 
with zero mass and no interactions with the outside world. Using dimensional reduction we 
obtained a two-dimensional equation for a massive spin-1/2 particle in curved space, which 
interacts via the gravitational metric g^ v and the electromagnetic potential A^. 

We also began with a trivial flat space metric in three-dimensions, Equation (J3J) and 
using dimensional reduction, derived a two-dimensional curved metric, Equation (1491) . and 
an electromagnetic vector potential, Equation fl50|) . 

Equations (1481) . ( )49l) and ( J50l) represent the main result of this paper. They provide a 
quantum mechanical treatment of a combined system of the Dirac Equation, electromag- 
netism and gravity, albeit in a low dimensional setting. We note however that this whole 
procedure can be carried out in higher dimensions. The only part of the derivation that is 
lacking in higher dimensions is an explicit solution of the Dirac spinor in terms of fractional 
derivatives like those given in Equations (1III Al) . If a higher dimensional version of this the- 
ory is constructed it would provide an alternate approach to current theories of quantum 
gravity. 



VII. CONCLUSIONS 



We have taken a model of an elastic medium and derived an equation of motion that 
has the same form as Dirac's equation in the presence of electromagnetism and gravity. We 
derived this equation by using the formalism of Cartan to reduce the quadratic form of 
the wave equation to the linear form of Dirac's equation. We showed that the dimensional 
reduction technique from Kaluza Klein theory produces not only the Einstein Field and 
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Maxwell's equations but also induces both mass and interaction terms into Dirac's Equation. 
The formalism demonstrates that a quantum mechanical treatment of the Einstein-Dirac- 
Maxwell equations can be derived from the equations of motion of the Fourier modes of an 
elastic solid and provides a new approach to theories of quantum gravity. 



APPENDIX A: QUANTIZATION COEFFICIENTS 

In diagonalizing the Hamiltonian given in Equation ( fl5l) we first write the field operators 
in terms of an intermediate set of ladder operators 

P = — (a ] -a 
n ' q \pl ^ 9 

— ( t 

Now define the vectors 

Q q = (ai,q, Cl2,q, ^3,5, CLl,-q, Cl2,-q, a 3,-q, <A,qj a 2,qi a 3,qi a l,-q) a \-q) a 3,-q 
Q\ = (°1,<?J a 2,<?5 a 3,q> a l,-g; a 2,-g; a 3,-q> a 2,c?) a 3,q> a l,-qi a 2-qi a 3,-qj 

This allows us to write H q = QlAQ q where the matrix A is given by 



.4 



T S 
S T 



16(A + 2/i) 
4(A + 2/i) 



and the nonzero elements of the matrices T and S are 
T 3 

Tx, 5 

^2,2 

^2,3 
^2,4 

^3,2 



4(A + 2^) 

8 (2q 2 2 + gf) /i 3 + 2/i + A (4 {Aql + qj) /i 2 + 1) 
16/i(A + 2/i) 

/i(3A + 2/i)g 2 g 3 
4(A + 2/i) 



4(A + 2/i) 
/i(3A + 2/i)g 2 g 3 
4(A + 2/i) 
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8 (qj + 2gf) /i 3 + 2/i + A (4 (ql + 4g 3 2 ) /i 2 + 1) 
16/i(A + 2//) 

4(A + 2/x) 
4(A + 2/x) 

__^3_ 

4(A + 2/i) 
3 

16(A + 2/i) 
4(A + 2/i) 

8 (2g 2 + gj) /i 3 + 2/i + A (4 (4g 2 + gf) /i 2 + 1) 
16/i(A + 2/t) 

/i(3A + 2/i)g 2 g 3 
4(A + 2/i) 

/*g3 

4(A + 2/i) 
/i(3A + 2/i)g 2 g 3 
4(A + 2/i) 

8 (g 2 + 2gj) /i 3 + 2/i + A (4 (gj + 4g 2 ) /i 2 + 1) 
16/i(A + 2/i) 

(2A + 2/i)g 2 

8(A + 2/i) 
(2A + 2/i)g 3 

8(A + 2/i) 
3 

~16(A + 2/i) 
(2A + 2/i)g 2 
8(A + 2/i) 

2/i (4/i 2 (2g 2 + gj) - 1) + A (4/i 2 (4g 2 + gj) - 1) 
16/i(A + 2/i) 

/i(3A + 2/i)g 2 g 3 

4(A + 2/i) 

(2A + 2/i)g 3 

8(A + 2/i) 
/i(3A + 2/i)g 2 g 3 

4(A + 2/i) 

2/i (4/i 2 (gj + 2gf) - 1) + A (4/i 2 (gj + 4gf) - 1) 
16/i(A + 2/i) 

3 

~16(A + 2/i) 
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5*4,11 

5*4,12 
55,8 
55,9 

55,10 
56,8 
56,9 

56,10 



(2A + 2/i)g 2 

8(A + 2/i) 
(2A + 2//)g 3 

8(A + 2/i) 

2/z (4/i 2 (2gj + gf ) - 1) + A (4/i 2 (4g 2 + g 3 2 ) - 1) 
16/i(A + 2/i) 

/i(3A + 2/i)g 2 g 3 

4(A + 2/i) 
(2A + 2/i)g 2 
8(A + 2/x) 
/i(3A + 2fx)q 2 q 3 

4(A + 2/i) 

2/i (4/i 2 (ql + 2g 2 ) - 1) + A (4/i 2 (gf + 4g 2 ) - 1) 



16/i(A + 2/i) 



(2A + 2/i)g 3 



8(A + 2/i) 

The methods outlined by Tikochinsky and Tsallis^ 1 ^ can now be applied to the matrix A. 
The final results allow us to define six creation operators b\ and six annihilation operators 
hi that satisfy 

Kq, b ],q>] = i$i,A,q' [Kq, b i,q'} = [ b h> b lq'} = 

These operators are Eigenstates of the Hamiltonian with eigenvalues 

1 



E 1>q 

E 2 ,q 



, v 12 + ql 

\\!<&+<ii 



E. 



3,g 



E 



Lq 



E, 



5,g 



E 



6,q 







//i(A + 2/i) 2 (2A + 2/i) 




(A + 2/i) 2 




2V2 X 


//i(A + 2/i) 2 (2A + 2/i) 




(A + 2/i) 2 




2V2 X 


//i(A + 2/^) 2 (2A + 2/i) 




(A + 2/i) 2 


\ 1+ 


2V2 X 


//i(A + 2/i) 2 (2A + 2/i) 




(A + 2/i) 2 



gf + ?l 



92 +<?! 



Q2 + Q3 



ql + ql 
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and the Hamiltonian has the diagonal form 



i 

The field displacement operators, Uij can be written in terms of these ladder operators. 
Denote the vector of field and ladder operators as 

X q = (Pl,q, Pl,qi Pj,,qi Ul,q, U 2 , q , U^ q , Pl-g, Pl-qi Pz-qi U l-qi U 2,-q, u 3,-q) 
B q = (b 1: q, b 2 ,q, 63,9, h,q, h,q, h,q, &l,<p &2,g> &3,g> &4,g> b\ q , ^ . 



This allows us to write 



Xi,q — ^ c i,jPi,q 



where the coefficients Cij are listed below. In writing the Cij coefficients we have defined the 
following quantities 



A 2 = 

A 3 = 

dx = 

d 2 = 

d 3 = 

d 5 = 

ni = 

n 2 = 



v 7 ^ vV(A + 2/i) 2 (2A + 2/i) 
V 



2V2~Jfi{\ + 2/i) 2 (2A + 2 fj) 



(A + 2/iy 



2y/2 J fi(\ + 2/i) 2 (2A + 2 fj) 



\] ' (A + 2/i) 2 +1 
8/i(2A + 2/i) (l6EqV + 16Eq 2 A/i 2 - 2/i - A) 

+ a (32Eq 2 /i 3 + 32Eq 2 A / u 2 - 10/i - 2A) 
8/i(2A + 2/i) (l6EqV + 16Eq 2 A / u 2 - 2/i - A) 

+ a (-32Eq 2 // 3 + 10/i + A (2 - 32Eq 2 /i 2 )) 
a (32Eq 2 /i 3 + 32Eq 2 A/i 2 - 10/i - 2A) 

- 8fi(2\ + 2fi) (l6Eq 2 /i 3 + 16Eq 2 A/i 2 - 2/i - A) 
12a/i + (A + 2/i) (64Eq 2 /i 4 + 96Eq 2 A/? + 4 (8Eq 2 A 2 + 3) fi 2 + I8\fi - 2A 2 ) 
(A + 2/i) (64Eq 2 // + 96Eq 2 A / u 3 + 4 (8Eq 2 A 2 + 3) fi 2 + 18A/i - 2A 2 ) - 12a/i 

2^2, 



g (gf + gf) 3/2 
^g 2 2 (2/i^-l)2 



2^ 



(gf + gf) 3/2 

^g 2 2 (2/i^-l) 2 
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n 3 = 8v / 2^8 / u 2 (2A + 2/i)^(2(2/i(A + 2 / u)A^ + 2/i(3A + 4 / u)^ 

+A 2 (16^V + 2/2 + A(16/Aj 2 + l))u q )a 2 - (A + 2/i)(4/i(A + 2/x)(2A + 2{i)\ 2 lu 2 
+4/i(2A + 2/i)(A + lO/i V! + (A + 2/i)A 2 (2/i(16/A; 2 + 5) + A(32/Aj 2 + 2))u q )a 
+4/i(A + 2/i) 2 (2A + 2 / u)(2/i(A + 2/i)A 2 .u; 2 + 2/x(3A + 4/i)u; 2 



n 4 = 8^- / u(A + 2/i) 2 (8(A + 2 / u)(2A + 2/i)(4AcjJ-3A 2 ^)/i 2 

+ a(A 2 c^((32/Ac> 2 + 2)A 2 + 6/i(16/A> 2 - 3)A + 4 / u 2 (16/Aj 2 - 3)) 

\ 1/2 / 

- 8(A - 6/j)h(2\ + 2fi)u 2 q ))j /UUafjL + (A + 2 /i) ((32 ^ - 2)A 2 

\ 1/2 

+ 6/i(16/A; 2 + 3)A + 4n 2 (16nW q + 3))) 2 



n 5 = 8v / 2^ 2 (2A + 2/i)cjJ(2(2 / u(A + 2/i)A 2 wJ + 2 / u(3A + 4 / u)^ 

+ A 3 (16^/i 3 + 2/i + A(16/AcrJ + l)K)a 2 + (A + 2/i)(4/i(A + 2//)(2A + 2/i)A 2 wJ 
+ 4/i(2A + 2/i) (A + lO//)^ 2 + (A + 2/i)A 3 (2/i(16/A; g 2 + 5) + \(32fi 2 u 2 + 2) V> 
+ 4/i(A + 2/i) 2 (2A + 2fj)(2fj,(X + 2/i)A 2 u; 2 + 2/i(3A + 4/i)cjJ + A 3 (16c^/i 3 + 2/i 

\l/2 / 

+ A(16/i 2 ^ 2 + 1)H))J /(JA + 2 / u)g 2 (8/i(2A + 2/i)(16u; 2 /i 3 - 2/i 

\ 1/2 

+ A(16/ArJ - 1)) + a(32^ g V - 10/i + A(32/i 2 ^ 2 - 2))) 2 



n 6 = 8^/i(A + 2/i) 2 (8(A + 2 / u)(2A + 2/i)(3A 3 cj g -4A^ 2 )/i 2 

+ a(A 3 u; g ((32/A; 2 + 2)A 2 + 6/i(16/i 2 u; 2 - 3)A + 4/i 2 (16/A; 2 - 3)) 

\l/2 / 

- 8(A - 6/i)/i(2A + 2fi)u 2 ))) /(((A + 2/i)((32/i 2 wJ - 2)A 2 + 6/i(16// 2 a; 2 + 3)A 

\ 1/2 

+ 4/i 2 (16/i 2 cuJ + 3)) - 12a/i) 2 



\l/2 / 




\ 1/2 




Finally, the nonzero components of the matrix c^- are 



Cl,3 = 



tv / 2cj g (a(2A-6^)(A+2^)A2-4^(2A+2^)(4A^(A+2^)-a(A-6M))^g) 
(A+2/i)d 2 re3<33 




iV2((A+2^)(-2A 2 +18MA+12^ 2 +4^((A+2^)(2A+2^)-2a)A 3 w 9 )-12o/i) 
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iV2g :) 

ty^(A+2»)g2((A+2»)(-2AA2+6»A2+4»(2A+2»)^ g )-8a»^ q ) 
tv / 2(A+2^)< ? 2(8a^ g +(A+2^)(-2AA 3 +6MA 3 +4/j(2A+2^)^ g )) 

d 5 HQU) q 



n 1 q 2 -2fj,n 1 q2UJ q 

tv / 2q2(a(2A+10/j+4/j(2A+2 At )A2^)-8M(A+2^)(2A+2/j)(2^A 2 <^ g +l)) 

d2«3<33 

ty / 292(8^(A+2^)(2A+2^)(2^A 3 a; g +l)+a(2A+10^+4^(2A+2 At )A3^)) 



n 2 — 2fj,n2Uq 

tv / 2(A+2/j) g 3((A+2/j)(-2AA2+6^A2+4^(2A+2^)^ g )-8a^q) 
iv / 2(A+2^)< ? 3(8a^ g +(A+2^)(-2AA 3 +6MA3+4/j(2A+2^)^ g )) 

ni—2finiuj q 

■iy / 2(a(2A+10^+4^(2A+2^)A2^ 9 )-8At(A+2^)(2A+2 At )(2 At A 2 ^+l)) 

ty / 2(8^(A+2^)(2A+2^)(2^A 3 a; g +l)+a(2A+10^+4^(2A+2^)A 3 a;q)) 

din 5 

4v / 2»(A+2u)^((A+2»)(2A+2»)(4»LJ 9 -A2)-2aA2) 
4 v / 2»(A+2u)^(2aA3 + (A+2»)(2A+2»)(4»^ g -A 3 )) 

4v / 2«(2A+2u)a;2 (4uA 2 w 9 a+a-4u(A+2^)) 
d2n 3 q 3 

4v / 2»(2A+2^)^(4 u A 3 w 9 a+a+4u(A+2u)) 
dl«593 

2\f2~iiqzw q 



niq2-2(j,niq2UJ q 

4 v / 2»(2A+2»)g2^(a(A 2 +4uc 1 ; g )-4»((A+2»)A2+4u(2A+2»)^)) 
^3^393 

4y / 2»(2A+2»)g2^(a(A 3 +4uc 1 ; g )+4»((A+2»)A 3 +4»(2A+2»)^)) 
din 5 q 3 

2y/2/j,q3Lii q 
n 2 q 2 — 2fj,n 2 q2UJ q 

4V2u< ?2 (2aA+(A+2»)(2A+2»)(-A+6H+4»(A+2u)A2^ g )) 
d4«4 

4-^/2^2 ((A+2u) (2A+2») ( - A+6u+4»(A+2u) A 3 a;q) -2aA) 
2v / 2"^g 



ni — 2[in\w q 

4y / 2»(2A+2»)^(a(A2+4^)-4»((A+2»)A2+4u(2A+2»)^ 9 )) 
(i 3 n 3 

4v / 2»(2A+2»)^ g (a(A 3 +4^ g )+4»((A+2»)A 3 +4u(2A+2»)^ g )) 

2\f2lLLO q 

n 2 — 2fm2U q 

4y / 2»< ?3 (2aA+(A+2 At ) (2A+2^) ( - A+6»+4»(A+2») X 2 ui q )) 

4y / 2»<? 3 ((A+2u) (2A+2») ( - A+6»+4»(A+2u) A 3 cj g ) -2aA) 
d 5 n6 
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iV2(l2a^+(A+2^)(-2A 2 +18MA+12^ 2 +4^(2a+(A+2/i)(2A+2/i))A 2 w g )) 

^4/24 

iV2((A+2^)(-2A 2 +18MA+12^ 2 +4^t((A+2^)(2A+2^)-2a)A 3 w g )-12a/i) 

ty / 2c 1 ;q(a(2A-6^)(A+2^)A2-4^(2A+2^)(4A^(A+2^)-a(A-6At))^) 
(A+2/i)d2«.3'23 

iy / 2^ g (a(2A-6^)(A+2^)A 3 +4/j(2A+2^)(a(A-6^)+4A^(A+2^))^ g ) 
(A+2/i)din5Q3 

iv2g3 

ty / 2g2(a(2A+10^+4^(2A+2^)A2^ g )-8M(A+2^)(2A+2^)(2^A 2 ^+l)) 

(fe«393 

ty / 292(8M(A+2 At )(2A+2^)(2^A 3 a; g +l)+a(2A+10^+4^(2A+2 At )A 3 ^)) 

di«5<?3 

i\/2g3 

ri2<?2-2/in 2 g2^9 

tv / 2(A+2/j)q2((A+2/j)(-2AA 2 +6/jA2+4^(2A+2 At )^ g )-8a^q) 
tv / 2(A+2^)< ? 2(8a^ 9 +(A+2^)(-2AA 3 +6MA3+4/j(2A+2^)^ 9 )) 

?V2 
ni— 2/j,niLO q 

tv / 2(a(2A+10^+4^(2A+2^)A2^ 9 )-8At(A+2^)(2A+2 At )(2 At A2^+l)) 

d2"3 

tv / 2(8^(A+2^)(2A+2 At )(2/jA 3 a; g +l)+a(2A+10^+4^(2A+2^)A3a; g )) 

din 5 

n2—2firi2UJq 

tv / 2(A+2/j)< ?3 ((A+2^)(-2AA 2 +6^A2+4^(2A+2^)^ g )^8a^q) 
(14714^ 

iv / 2(A+2/j)< ?3 (8a^ g +(A+2^)(-2AA 3 +6MA3+4/j(2A+2^)^ g )) 

d^nQLOg 

4v / 2M(2A+2^)w 2 (4/iA 2 w ? a+a-4^(A+2^)) 
d2"3<?3 

4 v / 2M(2A+2/i)a; 2 (4^A 3 w 9 a+a+4^(A+2^)) 

4y / 2^(A+2 A t)^ 9 ((A+2^)(2A+2/j)(4^a; g -A2)-2aA2) 
c(4n4 

4%/2^(A+2^)^(2aA3 + (A+2^)(2A+2^)(4^ g -A 3 )) 

2\/2^q 3 uj q 
rt2cj2-2/in2g2<^9 

4v / 2 W 2(2aA+(A+2^)(2A+2^)(^-A+6At+4^(A+2^)A 2 a;g)) 
C?4?24 

4 v / 2^g2((A+2/j)(2A+2^)(-A+6^+4^(A+2/j)A 3 a; g )-2aA) 

2y / 2/jg3i 1 ; g 
rii<j2-2/inig 2 ^9 

4v / 2/j(2A+2/j)<?2^ 9 (a(A2+4^ g )-4^((A+2^)A2+4^(2A+2/j)c 1 ; g )) 
rf3"393 

4v / 2^(2A+2^)< ? 2^(a(A3+4^ g )+4^((A+2^)A 3 +4^(2A+2^)c 1 ; g )) 
dm 5 q3 

ri2~2fin2U> q 

4v / 2^g3(2aA+(A+2/j)(2A+2/j)(-A+6At+4^(A+2^)A2^ g )) 
C?4?24 
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iV^iMqs ( ( A+2/i) (2 A+2/i) ( - A+6/i+4/i( A+2/i) A 3 w g )-2oA) 
C 12,6 - ^ 

_ 4v / 2^(2A+2/j)a; g (a(A2+4^ g )-4^((A+2^)A2+4^(2A+2/j)a; g )) 
12 > 9 d3n,3 

4y / 2^(2A+2/j)c 1 ; g (a(A3+4^ g )+4^((A+2^)A3+4^(2A+2^)a;q)) 
C 12,ll ~~ cijng 



APPENDIX B: AUXILIARY EQUATION 

We wish to show that 



dx 

We first consider the equation 



- + 7"r* - r„7" + 7"r„ = o 



BY 



dx v 

true in the unprimed coordinate system. But since the unprimed coordinate system is 
Euclidean space, the Christoffel symbols are identically zero. This allows us to write 

Since this is a tensor equation true in all frames, in the primed coordinate system we can 
immediately write 

Ki" + i^l = o 

Using 7 /M = S^S -1 , we have 

Wrs- 1 ) + (s^s^wl = o 

or 

(^s)rs- 1 + s&Yjs- 1 + sn^s- 1 ) + (s^s- 1 )^ = o. 

Multiplying by S^ 1 on the left and S on the right yields 

s-\& v s)f + + r&s- i )s + = o 

Finally, using Y u = (& v S~ l )S and again noting that d' u S~ l S = —S~ l d' v S we have, 

7*r„ - + (d' u r + 7%^) = (Bl) 
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APPENDIX C: DECOMPOSITION OF f 3 

In decomposing the matrix 7 /3 we seek a matrix 7j_ that satisfies 

{7 1 ,7±} = 0, {7 2 ,7±} = 0, 7l = l (CI) 
where from Equations (|4T|) and (!38|) we can write 

7 1 = ^((l-^ + ^f + ^h 3 )^ 



f = 5- 1 (u? 7 1 + (l-^)7 2 + ^7 3 )^ 



We can treat the matrices 7 M as vectors and use the cross product formula to compute an 
orthogonal vector. In other words we can write 



using the rules 



7± = 7 1 x 7 2 



7 1 x 7 2 = 7 3 7 2 x 7 3 = 7 1 7 1 x 7 3 = — 7 2 



This gives 

v± = S- 1 [\u\ul - u\{\ - ul)}^ + [u\u\ - (1 - u\)u 2 ^ 2 

+ [{l-u\){l-ul)~u l 2 u\] 1 ^S 

with 

v± 

7i = I — 1 
\v±\ 

and \v± \ = Jvj_. It can be verified directly that 7j_ satisfies Equation fIClj) . 



APPENDIX D: FOURIER COMPONENTS 

We wish to relate the Fourier components of the field displacements Ui, when expanded 
in the basis e iq ' x ' to those expanded in the basis e tqx . The two coordinate systems are related 
by u^ = x^ — x'^ where the field displacements Ui are assumed to be small. We write 



iqx 



iq-(x'+u) 
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~ ^2u m e igx '(l + iq-u) 
i 

q k 

» Ev ,?1 '( 1 +«E?' if / 1 ') 

<? k 

= Y. u m( 5 kfi + n-u k )e l{q+k) - x ' 

qk 

correct to second order in u^ q . 
Setting q + k = q' we have 



m m = H^W-fc^M + l (q' ~ k ) ■ Uk)e ig '" x ' 

q'k 

so that the components of the field displacements in the primed frame are 



u m ' = u m'-k{kfi + %(rri - k)u k ) 
k 

We can now relate the energy eigenstates which are expressed in the unprimed frame to the 
Fourier components in the primed frame. In a low energy theory in which only uq is present, 
to second order in Ui, the only Fourier component in the primed frame that is nonzero is 
u m i =0 . If there is enough energy to excite the fields u and u±l. Then the only nonzero 
modes in the primed frame are u m > =0 and u m '=±i and so on. 
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